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1 Introduction 



A future linear collider will be an excellent tool to search for and investigate extensions 
of the Standard Model (SM) of particle physics Q- One particularly attractive extension of 
the SM is Supersymmetry (SUSY) 0, which solves several conceptual problems of the SM. 
Apart from their direct production, also virtual effects of SUSY particles may lead to observable 
deviations from the SM expectations. In particular, top quark pair production at a linear collider 
may be a sensitive probe of such effects. Very high energy scales are involved in the production 
and decay of top quarks. Moreover, since they decay very quickly, the spin of top quarks is 
not affected by hadronization effects and becomes an additional observable to probe top quark 
interactions. At a future linear e + e~ collider, the electron (and possibly also the positron) beam 
may have a substantial longitudinal polarization, which will be an asset to study top quark spin 
phenomena. We therefore study in this paper the impact of virtual effects of SUSY particles 
on spin properties of tt pairs in collisions. We restrict ourselves here to the SUSY QCD 
sector of the Minimal Super symmetric Standard Model (MSSM). SUSY QCD corrections to the 
(spin-summed) differential cross section for — > tt have already been studied quite some 
time ago 0, and we extend these results by keeping the full information on the tt spin state. 
The full MSSM corrections to the spin-summed differential cross section have been calculated 
in a. 

In section 2 we define the spin observables that we calculate in this paper and also discuss how 
they can be measured. Section 3 gives analytic results for these observables, and section 4 
contains numerical results for specific choices of the SUSY QCD parameters. In section 5 we 
present our conclusions. 



2 Spin observables 

We consider the reaction 

e+(p + ^ + ) + e -(p-,X-)^(y\Z*)^t(k t )+t(k)+X, (1) 

where A,_ (k + ) denotes the longitudinal polarization of the electron (positron) beam 1 . Within 
the Standard Model, spin effects of top quarks in reaction © have been analysed first in ref. 
0. QCD corrections to the production of top quark pairs, including the full information about 
their spins, can be found in refs. ABED- Fully analytic results for the top quark polarization |0 
and a specific spin correlation @ to order oc^ are also available. 

The top quark polarization is defined as two times the expectation value of the top quark spin 
operator S t . The operator S f acts on the tensor product of the t and t spin spaces and is given by 
St = ? <8> 11, where the first (second) factor in the tensor product refers to the t (f) spin space. 
(The spin operator of the top antiquark is defined by Sf = 11® S.) The expectation value is 

'For a right-handed electron (positron), X T = +1. 
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taken with respect to the spin degrees of freedom of the it sample described by a spin density 
matrix R, i.e. 

P, = 2<Sr>=2^l (2) 

For details on the definition and computation of R, see e.g. 0. The polarization of the top 
antiquark Pf is defined by replacing S t by S? in ©. For top quark pairs produced by CP invariant 
interactions, we have = P r . The spin correlations between t and t can be calculated by using 
the matrix 

Tr \RStiStj] 

Using arbitrary spin quantization axes a and b for the t and t spins, the spin correlation with 
respect to these axes is given by 

t «-AV(P,-8)(Prt) (4) 



l-(P,-8)Yl-(P<-«>) 2 

The directions a, b can be chosen arbitrarily. Different choices will yield different values for the 
spin correlation c(a,b). The spin properties of the top quarks and antiquarks can be measured 
by analysing the angular distributions of the t and t decay products. For example, if both t 
and t decay semileptonically, t — > bt+Vi, t — > bi'~\g, the following double differential lepton 
angular distribution is sensitive to the /J spin state: 

1 J 2 o 1 , 

= -(1 +Bicos0 + +B 2 cos0_ -Ccos0,cos0_) , (5) 

a<icos0 + <icos0_ 4 

with a being the cross section for the channel under consideration. In Eq. © + (0_) denotes 
the angle between the direction of flight of the lepton £ + (£'~) in the t (f) rest frame and the cho- 
sen spin quantization axis a (b). The coefficients B\ 2 and C are related to the mean (averaged 
over the scattering angle) t (f) polarization and spin correlation projected onto the directions a 
and b. Using the double pole approximation [ 10 1 for the t and t propagators, one obtains for the 
so-called factorizable contributions [fTTHT2l 



Bi = K + P f a, 



B 2 = — K_Pf • b, 



C = K+K-aiQjbj, (6) 
where the overline indicates the average over the scattering angle, e.g. 

J^dyTrlRSfk] 



P r a = 2 — ' / ' \ (7) 

etc., where y is the cosine of the top quark scattering angle. In ©, K± is the spin analysing 
power of the charged lepton £ ± . At leading order, K± = +1. QCD corrections to this result are 
at the per mill level [ 13] . SUSY QCD corrections to the spin analysing power K± are exactly 
zero [fill. 
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3 Analytic results 



We now turn towards the calculation of the SUSY QCD corrections to the polarization and 
spin correlations of top quark pairs produced in e + e~ collisions. These corrections directly 
determine the SUSY QCD corrections to the double lepton distribution © within the double 
pole approximation, since the corrections to the LO result K± = +1 are exactly zero and the 
non-factorizable contributions due to SUSY particles also vanish within that approximation. 

The amplitude for reaction (HJ) including SUSY QCD corrections may be written as follows: 

iT/i = ^|xWv(p+)(g^-^Y5)"(p-)^-v(p+)Y^(P-W^ (8) 

where g e v = — j + 2sin 2 iV, and g e A = —\, with i3w denoting the weak mixing angle. The 
function % is given by 

xOO = — Ti S ~' (9) 

4 sin f>vi/cos 2 f>iv s — m z 
where mz stands for the mass of the Z boson. We neglect the Z width, since we work at lowest 
order in the electroweak coupling and the cm. energy is far above mz- The hadronic currents 
have a formf actor decomposition as follows: 



H^ y =u(k t ) 



Vzjf-A z ,jfl5+Szx 



2m t 



v(k) (10) 



with 



Vz,y = Vl y + Vl v 

AZ,J = A Z y+A Z y. (11) 

In (fTTTl . Vy = Q t , where Q t denotes the electric charge of the top quark in units of e = ^/Ana, 
Ay = 0, and V® = g' v = j — | sin 2 f>iy , A z = g' A = j are the tree level vector- and the axial- vector 
couplings of the top quark to the Z boson. 

The one-loop SUSY QCD contributions to the different form factors are denoted by Vy z , Ay Z 
and Sy,z- Scalar and pseudoscalar couplings proportional to (k t + kjY and (k t + kfYy^ have been 
neglected in ®, since they induce contributions proportional to the electron mass. In addition 
CP violating formfactors proportional to (k t — kjYj5 are possible in SUSY QCD through a 
complex phase in the squark mass matrices. In BT31 it has been shown that the dependence 
of the cross section on these phases is weak and that CP odd asymmetries are typically of the 
order of 10 -3 . We therefore set these phases to zero in the following. To make this paper self- 
contained we list the form factors Vy Z , Ay Z and Sy.z in the appendix. We have performed an 
analytic comparison to the corresponding results in and found complete agreement. 

We define 

/ll(lr) = -Qt+X(8v+8A)(8 t v ± 8 t A)^ 

/rr(rl) = -Qt+X(g e v-g e A )(gvTg A ), (12) 



3 



and 

P± = 1-X-X+±(X--X+). (13) 
The electroweak couplings that enter the Born results are then given by 

8VV = l[P+(fRR + fRL) 2 ±P-(fLL + fLR) 2 }, 

S^A = l[P+(fRR-fRL) 2 ±P-(fLL-fLR) 2 }, 

SVA = l[P+(fRR-fRL)±P-(flL-fl R )]- d4) 

Likewise, defining 



gLL(LR) = l(g e y+g e A )(V^±A l z )-(V^±A\) 

8rr(rl) = l(g e v -g e v )(V^TA l z )-(V^TA\) : 

s L (R) = x(g e v ±g e A )S z -Sy,. (15) 

the SUSY QCD contributions may be written in terms of the following quantities: 



hyy = £ [P+ {fRR + fRL) (gRR + gRL) ± P- (Al + fut) (gLL + glJl)] , 
h AA = ^[P+(fRR-fRL)(gRR-gRL)±P-(fLL-fLR)(gLL-gLR)}, 

Re hy A = ^Re [P + (f RR g RR - f RL g R L) ± P- (/ll8ll - /zjegzje) ] , 

Im h y A = ^Im [P + (f RL g R R - fRRgRL) ± P (flRgLL ~ fLLgLR)] , 
4 = \[P+(fRR + fRL)sR±P-(fLL + fLR)s L }, 

4 = -\[P + {fRR-fRL)SR±P-{fLL-fLR)SL]. (16) 

It is convenient to write the results in terms of the electron and top quark directions p and k 
defined in the cm. system, the cosine of the scattering angle y = p • k, the scaled top quark mass 
r = 2m t I y/s, and the top quark velocity (3 = \J\-r 1 . 

The differential cross section including the SUSY QCD corrections reads: 

+ p 2 (l+/)(^+ 4 + 2Re/i+ 4 )+4pv(^+ A + 2Re^-2p 2 (l-/)Re4|, (17) 



where 



4na 2 



(18) 



and dc°/dy is obtained by setting hy V = hj^ = hy A = Sy = 0. We further introduce a vector 
perpendicular to k in the production plane = p — yk and a vector normal to this plane, 
n = p x k. The top quark polarization including the SUSY QCD corrections is equal to the top 
antiquark polarization and reads: 

P, = Ff + P} 

3Af c p 



> P t- 



+ r 



W +y 2 ) {Sva + 2Re Ka) +y (Svv + 2R e Kv) + tfy (Saa + 2Re h A^)] fe 

P 2 



fry {g VA + 2Re Ka) + Sw + 2Re h vv - T ( Re s v - $y Re S A ) 



+ 



B 2 

2pYlm h+ A + Z- (yhn s+ - Blm s+) 



n 



dy 



-l 



)do l ( do 



dy \ dy 



-l 



(19) 



For the matrix Qj defined in © we find 

^0 ,A 1 5 [ d(5 l fdO°Y r 



+ c 



+ 



^ {%) 1 { +2Re " p2 ^ +2Re [*W " I 8 

(y 2 + B 2 (l - y 2 )) (g+ v + 2Re h+ v ) + B 2 y 2 + 2Re *+J + 2By (#+ A + 2Re *+J 



+ 2B 2 (l-y 2 )Res 



+ 



A A 1 ^ 



+ r 



P 2 



y (g+ v + 2Re + B (g+ A + 2Re h+ A ) - % (jRe 4 - BRe s +) 



A^ - Ay ~\~ Ay A^' 



+ 2Blm/z 



VA 



k t rij + k; rij 



2yrlm h yA + - (im s v - Bylm s A ) 



da 1 / Jo 
\~dy 



(20) 



The Born results Pf and Cj^ are obtained from (fT9l and d20l by setting /z^ A = = = s 
s A = do 1 /dy = 0. 



For fully polarized electrons (or positrons) a so-called 'optimal spin basis' can be constructed. 
This is an axis d with respect to which the t and t spins are 100% correlated at the tree level in 
the Standard Model for any velocity and scattering angle [ 16 1. This axis d is the solution of the 
equation 



(21) 



One gets 

d = xk+ v / l-x 2 k ± , (22) 

with x E [—1,1] only if either P + = or P = 0. For P + = 0, which can be realized with 
left-handed electrons (k- = — 1), one finds 

x= fLL($+y)+fLR(y-V) 

[(1 +yV) 2 fh+ (1 -^)% + 2(j2p2 + ! _2 P 2 )/ll/ ^] 1/2- 

For right-handed electrons, the optimal basis is obtained by replacing fn — > /rr, /lr — ► //jz, 

. . -p^o 

in Eq. (I231) . Note that at threshold d — >p, i.e. the optimal basis at threshold is defined by the 

direction of the beam, while in the high-energy limit d — >k, i.e. the optimal basis coincides 
with the helicity basis. By analytically evaluating d[C\dj we find that the virtual SUSY QCD 
corrections to the it spin correlations in the optimal basis are exactly zero. 



4 Numerical results 

In this section we present numerical results for the SUSY QCD corrections to the top quark 
polarization and tt spin correlations. We also include a discussion of the corrections to the 
differential cross section and compare our results to the literature. 

We take into account the effects of mixing of the chiral components of the top squark. The stop 
mass matrix can be expressed in terms of MSSM parameters as follows: 



M| + m 2 + m|(i-a4)cos2p m t (A t - p cot $) \ 

m,(A,-jUCOtp) M 2 7 +m 2 + m|2^ 2 y cos2p J ' U> 




where Mq, Mq are the soft SUSY-breaking parameters for the squark doublet qi (q = t,b) 
and the top squark singlet Ir, respectively. Further, A t is the stop soft SUSY-breaking trilinear 
coupling, and fj is the SUSY-preserving bilinear Higgs coupling. The ratio of the two Higgs vac- 
uum expectation values is given by tanp\ and we use the abbreviation sw = sinG^. The squared 
physical masses of the stops are the eigenvalues of the above matrix. In order to simplify the 
discussion, we set tan (3 = 1 for all following results. Further, we assume that the sbottom mass 
matrix is diagonal with degenerate mass eigenvalues, fAf- 2 = diag(m~,m~). Neglecting mt, in 
the sbottom mass matrix this leads to Mq = my, and the stop mass matrix simplifies under the 
above assumptions to 

= { m t M LR Ml+mf J' (25) 

with Mlr =Af—fi. The stop mass eigenstates are obtained from the chiral states by a rotation: 

h \ = f cos0 ? sin0 f \ / t L \ 
h J V -sin0 f - cose,- )\h )' 
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M LR [GeV] 

Figure 1: Relative correction to the total cross section as a function of the mixing parameter 
Mlr for a fixed mixing angle 0^ = rc/4. Shown is the correction for two different gluino masses 
m~ g = 150 GeV and m, = 250 GeV. 



Maximal mixing (0f = | and M/j? 7^ 0) corresponds to Mq = mi. The latter relation will also 
be assumed for Mlr = 0, leading to the following stop mass eigenvalues 2 : 

m fl2 = ^mj + mj±m t M LR . (27) 

Note that we use here the same set of assumptions on the squark mass matrices as we did in our 
study of the SUSY QCD corrections in the decay of polarized top quarks ifbH . Further we use 
sin Qw = 0.2236, a s = 0.1 1, and we set the top mass to m t = 174 GeV and the sbottom mass 
that enters Eq. (E3> to m h = 100 GeV. 

Fig. 1 shows the relative SUSY QCD correction a 1 /o° to the total cross section for e + e~ — > tt 
with unpolarized beams at ^Js = 500 GeV as a function of the mixing parameter Mlr, where 
G° and a 1 are obtained from Eq. (fTTl) by integrating over y. Shown are the relative corrections 
for two different gluino masses, namely nig = 150 GeV and nig = 250 GeV. For a large mixing 
parameter Mlr and a small gluino mass of nig = 150 GeV we find a large negative correction. 
The correction decreases as the gluino mass increases. A mixing parameter of Mlr = 200 GeV 
corresponds to a light stop mass of mf 2 = 74 GeV, which is above the current experimental 
lower limit ifTTll . With our choice of the masses, we are far away from the threshold singularity 
at m t = nig + m~ u where a more sophisticated calculation is necessary. 

Fig. 2 shows the differential cross section do/dy, again for two different gluino masses nig — 
150 GeV and nig = 250 GeV, and for the cases of 'no mixing' (Mlr = 0) and 'mixing' (Mlr = 
200 GeV and 0^ = 7t/4), again at y/s = 500 GeV. We have compared our results for a and 
do/dy with and found agreement with their Fig. 3 (no mixing case), while we disagree with 
the results depicted in Fig. 4 (a and forward-backward asymmetry with stop mixing). We have 
also compared our results including the mixing with [4, 18 1 and find complete agreement. 

2 Note that by fixing f - = | the light stop can be either ti or I2 depending on the sign of Mlr. 
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-1 -0.5 0.5 1 

y 

Figure 2: Relative correction to the differential cross section dc/dy at y/s = 500 GeV for the 
cases of no mixing {Mm = 0) and mixing (0^ = 7l/4 and M LR = 200 GeV). 



We now turn towards the discussion of the SUSY QCD corrections to the tt spin properties. 

In Fig. 3 we investigate the expectation value of the top spin operator as a function of the 
centre-of-mass energy. We have computed the average projected polarization defined in Eq. Q) 
for three choices of the quantization axis a, namely for a = k (flight direction of the top), for 
a = p (electron beam direction), and for a = n (normal to the event plane). These quantities are 
shown in three different plots, where thin curves correspond to the tree level results and the thick 
curves are the relative corrections in percent. The corrections are shown for the case of mixing 
(Qq = 7t/4 and Mir = 200 GeV) and a gluino mass of trig = 150 GeV. For the polarizations of the 
initial beams we choose A,+ = and consider the three cases X_ = — 1 , 0, +1 . The projection of 
the top quark polarization onto n vanishes at tree level, and thus we only show the contribution 
from SUSY QCD absorptive parts in percent. In all cases SUSY QCD effects change the tree 
level results by less than 1% and vanish at threshold. 

In Fig. 4 we show the averaged spin correlations djCjjbj for the choices a = b = k (helicity cor- 
relation), a = b = p (beamline correlation), and a = p, b = k, for the same choice of parameters 
as in Fig. 3. Again the SUSY QCD correction are tiny. Fig. 5 shows the correlations for the 
choices a = k, b = n and a = p, b = n. The first of these two choices of spin quantization axes 
leads to SUSY QCD effects slightly larger than 1% around cm. energies of 700 GeV and for a 
fully polarized electron beam. 



5 Conclusions 

In this paper we have derived analytic expressions for the SUSY QCD corrections to the po- 
larization and spin correlations of tt pairs produced in e + e~ annihilation with longitudinally 
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[GeVl 




y/t [GeV] 

Figure 3: Average projected top quark polarization P r a defined in © for the choices a = k 
(top), a = p (middle), a = n (bottom) as a function of the centre-of-mass energy. In each plot 
we show the tree level results (thin lines) and the relative corrections in percent (thick lines) for 
unpolarized positrons and the three cases X_ = —1,0, +1. 
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Figure 4: Same as Fig. 3, but for the quantities ktCijkj (top), PiCtjpj (middle), and piCij 
(bottom). 
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Figure 5: Same as Fig. 3, but for the quantities kiCijhj (top) and PiCijhj (bottom). 
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polarized beams. The results depend in particular on the gluino mass and the masses of the 
scalar partners of the top quark. The latter masses depend on the mixing in the stop sector. For 
maximal mixing, the SUSY QCD corrections to the cross section are negative and reach values 
of about -1.3% (-5%) for a gluino mass of 250 GeV (150 GeV) and a light stop mass of 74 
GeV. For the same choice of parameters, the it spin observables typically receive SUSY QCD 
corrections well below 1%. 
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Appendix 

Here we list explicit results for the formfactors defined in Eqs. (fTTjb . (fTTT) . Apart from the Stan- 
dard Model parameters m t , a s , and electroweak couplings defined in section 3, the formfactors 
depend on the gluino mass mg, the masses of the two physical top squarks , and the mixing 
angle df that determines how the top squark mass eigenstates are related to the weak eigenstates, 
cf. Eq. 



We find (with C F = (N% - l)/2/N c = 4/3): 



V, 1 = ^C F Q t [C\\ +C%] +Q t ^±^, (A.1) 



V} = —C F [ (g A cos 2 f - Q t sin 2 tV) C l 2 \ + (g' A sin 2 8? - Q t sin 2 C 2 2 |] 

{ 8Z^ + 8Z L t 8Z R - 8Z L 
+ Sv o Sa o > ( A - 2 ) 



A] = ^C F Q t [C\\-C%\ cos29 f -& 5Z * 2 5ZL , (A.3) 

A\ = ^C F |2[(^cos 2 e r -e f sin 2 T3^)cii-(^sin 2 f -asin 2 T3,y)c| 2 ]cos20 f 
+ ^sin 2 20KC- + C 2 i)}-^^ Z ^ + ^? Z ^, (A.4) 
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Sy = —C F Qtm t 



Sy ~\~ Sy 



(A.5) 



where 



11(22) 

Jy 



m t (c;; (22) +C^ 1(22) ) ±m |S in2e f - (c^ +€[1^) , (A.6) 
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S z = ||c F m f 1 2 cos 2 r -ft sin 2 r> M /)4 1 + 2 sin 2 9 f - ft sin 2 r>iv)^ 

- g\ 8^20,-008 20^ (cj 2 + C} 2 + Cf +Cf{) J. (A.7) 

In the above expressions, the one-loop integrals Cq , . . .C 2 4 are defined by the decomposition of 
Passarino and Veltman ifT^I . 



^ nt 2 7 r/ 2 — mi + fe] [(/ — A, ) 2 — + ie] [(/ + *r) 2 — mi + fe] 



with {ky = k,+kt): 



Cj — ^f/jCy + ky/jC^i 

C l /lv = kt/jktvC 1 ^ + ky p ky v C~/ 2 - (k t/J k Vv + k tv k VlJ ) C 2 J 3 + g^C{ A . (A.9) 

The quantities 8Z#x denote the one-loop renormalization constants for the chiral components 
of the top quark field in the on-shell renormalization scheme. They are given explicitly by 

&m = ^{2mj[(B\y+(Bly]+2m § m t sm2Q f [(Bl ) y~(Bly] 



+ 5}+5 2 ±cos20 f [B\-B\] }, (A.10) 



where 



( ^) 4 - d f ddl Ml 

in 2 J [/ 2 -m| + z'e][(/-^) 2 -m? + z'e] 



4„= HSr- " M .j., M ,„'',« j., M =K;(-k„)B\, (A.ii) 



and 



(^1))' = ^ W (A. 12) 
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